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A total of 614 upward through-going muons of minimum energy 1.6 GeV are observed 
by Super-Kamiokande during 537 detector live days. The measured muon flux is (1.74 ± 
0.07(stat.) ± 0.02(sys.)) x 10~ 13 cm~ 2 s _1 sr _1 compared to an expected flux of (1.97 ± 
0.44(theo.)) x 10 _13 cm — 2 s _1 sr — The absolute measured flux is in agreement with the 
prediction within the errors. However, the zenith angle dependence of the observed 
upward through-going muon flux does not agree with no-oscillation predictions. The 
observed distortion in shape is consistent with the <-> v T oscillation hypothesis with 
sin 2 20 > 0.4 and 1 x 10~ 3 < Am 2 < 1 x lO^ 1 eV 2 at 90 % confidence level. 
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Energetic atmospheric v ^ or 9^ passing through the Earth interact with the rock surrounding the 
Super-Kamiokande ( "Super-K" ) detector and produce muons via weak interactions. While those neutrino- 
induced muons traveling downwards are impossible to differentiate from the constant rain of cosmic ray 
muons, upward-going muons are mostly or induced, because upward-going cosmic ray muons cannot 
penetrate through the whole Earth and v e and D e induced electrons and positrons shower and die out 
in the rock before reaching the detector. Those muons energetic enough to cross the entire detector are 
defined as "upward through-going muons" . The mean energy of their parent neutrinos is approximately 
100 GeV. Neutrinos arriving vertically travel roughly 13,000 km, while those coming from near the horizon 
originate only ~500 km away. 

Previously published results on atmospheric neutrinos with average energies below ~10 GeV have 
indicated an anomalously low v^/v e ratio Jl)-^ and have also reported a strong zenith angle dependence 
This has been interpreted as a possible signature of neutrino oscillations. Recent results from this 
experiment (5|(| have shown strong evidence for <-> v T oscillations . These results have reported on 
lower energy and v e neutrinos which interacted in the water of the detector itself, hereafter referred 
to as "contained" events. 

The oscillation hypothesis has also been suggested to explain the anomalous upward through-going 
muon zenith angle distributions observed by Kamiokande j|] an d MACRO || as well as the low absolute 
upward-going muon flux seen in MACRO. However, the absolute upward-going muon fluxes measured in 
Kamiokande, 1MB fTcf] , and Baksan were consistent with the no-oscillation expectations within the 
large errors present in the absolute flux predictions. 

We make the first report on the measurement of upward through-going muon flux and its zenith- 
angle distribution as observed by Super-K. The experimental site is located at the Kamioka Observatory, 
Institute for Cosmic Ray Research, the University of Tokyo, 1000 m underground in the Kamioka mine, 
Gifu prefecture, Japan. 

The Super-K detector is a 50 kton cylindrical water Cherenkov calorimeter. The detector is divided by 
an optical barrier instrumented with photomultiplier tubes ( "PMT" s) into a cylindrical primary detector 
region (the Inner Detector, or "ID") and a surrounding shell of water (the Outer Detector, or "OD") 
serving as a cosmic ray veto counter. Details of the detector can be found in reference ||. 

The cosmic ray muon rate at Super-K is 2.2 Hz. The trigger efficiency for a muon entering the detector 
with momentum more than 200 MeV/c is ^100% for all zenith angles. The nominal detector effective 
area for upward through- going muons with a track length > 7m in the ID is ~1200 m 2 . 

The data used in this analysis were taken from Apr. 1996 to Jan. 1998, corresponding to 537 days of 
detector livetime. Event reconstruction is made by means of the charge and timing information recorded 
by each hit PMT. The direction of a muon track is first reconstructed by several automated grid search 
methods, which find the track by minimizing the width of the residual distribution of the photon time- 
of- flight subtracted ID PMT times. Details of one such muon fitter are described elsewhere fL2"f . 

A minimum track length cut of 7m (~1.6 GeV) was applied. To reduce the abundant downward-going 
cosmic ray muons, events satisfying cos 6 < 0.1 are selected, where O is the zenith angle of the muon 
track, with cos O < corresponding to upward-going events. Muons which leave both entrance and exit 
signal clusters in the OD are regarded as through-going. After a visual scan by two independent groups 
(event loss probability < 0.01%) and a final direction hand-fit, 614 upward through-going muon events 
with cos 6 < are observed. Different hand fits are consistent with each other within 1.5°. They are 
shown to be unbiased in blind tests using Monte-Carlo (MC) simulated events, with deviations between 
the reconstructed track direction and the real muon direction (A8 rec ) estimated to be 1.4°. Using this 
same MC, the directional correlation between a muon and its parent neutrino is estimated to be 4.1°, 
including contributions from the muon production angle and from multiple Coulomb scatterings in the 
rock. 

Due to the finite fitter resolution and multiple Coulomb scattering in the nearby rock, some down-going 
cosmic ray muons may appear to have cos O < 0. Figure [l] illustrates the estimation of this contamination. 
Assuming this background continues to fall exponentially as cos O decreases, the contribution to apparent 
upward-going muons is estimated to be 4.3±0.4 events, all contained in the —0.1 < cosO < zenith angle 
bin. The contamination at the Kamioka site due to cosmic ray photoproduced upward-going pions [Of 
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meeting the 7m track length requirement is estimated to be < 1%. 

The total detection efficiency of the complete data reduction process for upward through-going muons 
is estimated by a Monte Carlo simulation to be >99% which is almost isotropic for —1 < cos 9 < 0. 
Using the upward/downward symmetry of the detector configuration, the validity of this Monte Carlo 
program has been checked by real cosmic ray downward through-going muons. 

This analysis used a model which is a combination of the Bartol atmospheric neutrino flux model 



and a neutrino interaction model composed of quasi-elastic scattering ]15| + single-pion production |16| 
deep inelastic scattering (DIS) multi-pion production. The DIS cross-section is based on the parton distri- 
bution functions (PDF) of GRV94DIS JlT]] with the additional kinematic constraint of W > 1.4 GeV/c 2 . 
Lohmann's muon energy loss formula in standard rock fl8| is then employed to analytically calculate 
the expected muon flux at the detector. This flux is compared to three other analytic calculations to 
estimate the mo del- dependent uncertainties of the expected muon flux. The other flux calculations use 
the various pairs of the Bartol flux, the GRV94DIS PDF, the atmospheric neutrino flux model calculated 
by Honda et al @, and the CTEQ3M [|o|]PDF. These comparisons yield ±10% of uncertainty for the 
absolute flux normalization and —3.7% to +1.6% for the bin-by-bin shape difference in the zenith-angle 
distribution. The shape difference is due mostly to the input flux models. 

The Bartol+GRV94DIS calculation results in an expected muon flux $ t heo of (1.97 ± 0.44(theo.)) x 
10~ 13 cm~ 2 s _1 sr _1 (cos 6 < 0), where the estimated theoretical uncertainties are described in Table [j]. 
The dominant error comes from the absolute normalization uncertainty in the neutrino flux, which is 
estimated to be approximately ±20% 0,|l9|,^l) for neutrino energies above several GeV. 

Given the detector live time T, the effective area for upward through-going muons 5(0), and the 
detection efficiency e(9), the upward through-going muon flux is calculated by the formula: 

« £ 1 11 

$obs = 2^ 



-£(6,) 5(9,) 2tt 
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where the suffix j represents each event number, 2ir is the total solid angle covered by the detector 
for upward through-going muons, N corresponds to the total number of observed muon events (614). 
Subsequently, we subtract the cosmic ray muon contamination (4.3 events) from the most horizontal bin 
(-0.1<cos9<0). 

Conceivable experimental systematic errors are summarized in Table [n| Including these experimen- 
tal systematic errors, the observed upward through-going muon flux is: & b s = (1.74 ± 0.07(stat.) ± 
0.02(sys.)) x lO^cm^s^sr- 1 . 

Fig. H shows the flux as a function of zenith angle. The shape of the distribution is not well represented 
by the theoretical prediction without neutrino oscillation, having a x 2 /degrees of freedom (dof) = 18.7/9 
corresponding to 2.8% probability. This shape comparison is done after multiplying the expected flux by 
a free-running normalization factor (1+a^), whose best fit value is a M = —14%. 

A set of neutrino oscillation hypotheses are then tested using the zenith angle distribution. The 
expected flux ( (cZ<& /<iSl) osc ) for a given set of Am 2 and sin 2 29 is calculated and the same binning 
(dcos9=0.1) is applied to this flux as to the data. To test the validity of a given oscillation hypoth- 
esis, we minimize a \ 2 which is defined as: 




where a s tat,i (&sys,i) is the statistical (experimental systematic) error in the observed flux (d$> / d£i) l obs 
for the iih bin, and (1 + a^) is an absolute normalization factor of the expected flux. The absolute 
flux normalization error a a ^ is estimated to be ±22 % by adding in quadrature the bin-to-bin correlated 
experimental errors and theoretical uncertainties in Table |. Based on the bin- by-bin correlated systematic 
errors in Table [n] added in quadrature, we estimate <y S y s ,i to range from ±(0.3—3.8)%. Then, the minimum 
X 2 (Xmin) is searched for on the Am 2 — sin 2 28 plane. 
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Assuming <-> v T oscillations, Xmin( = 7.5/8 dof) occurs at (sin 2 29, Am 2 ) = (0.95, 5.9 x 10 _3 eV 2 ) and 
= +12%, in good agreement with the overall normalization found in the contained event analysis Q, 
although the of this analysis refers to the flux normalization of neutrino energies predominantly around 
100 GeV. For the null oscillation case (sin 2 29=0), we obtain x 2 of 19.2 at a best-fit = —14% using 
the same x 2 definition. The zenith angle distribution of (1 + a il )(d^ / d£t,) l osc for the best fit parameters 
is shown in Fig. |^ together with the data. Figure || shows the confidence intevals on the (sin 2 29, Am 2 ) 
plane for <-> v T oscillations. The 90% C.L. contour marks the line of Xmin + 4-6- If we replace the 
Bartol neutrino flux Jl4j by Honda's jl9| and/or the GRV94DIS parton distribution functions |l7j by 
CTEQ3M p0| , the allowed region contours are similar to those presented in Fig. |3[ Consequently, we 
find that the zenith angle dependence is in favor of the «-> v T oscillation hypothesis and supports the 
Super-K contained event analysis . It is also consistent with the data presented in the Kamiokande § 
and MACRO |J upward-going muon analyses. Interactions of v T in the rock below is estimated at less 
than a few percent and neglected in this analysis. Oscillation of to v e in this range of parameter space 
has been ruled out by the CHOOZ experiment [^2| . 

In conclusion, based on 614 upward through-going muon events during 537 detector live days, the 
flux of the upward through-going muons (>1.6 GeV) is measured with the Super-K detector: $ t, s = 
(1.74±0.07(stat.)±0.02(sys.)) x 10 _13 cm~ 2 s _1 sr _1 . This is compared with the expected flux of <& t heo = 
(1.97 ± 0.44(theo.)) x 10 -13 cm~ 2 s -1 sr -1 . The absolute observed upward through-going muon flux is in 
agreement with the expected flux within the relatively large uncertainties in the theoretical calculations. 
We find that the zenith angle dependence does not agree with the theoretical expectation without neutrino 
oscillations at the 97% C.L. However, the «-> v T oscillation hypothesis with sin 2 29 > 0.4 and 1 x 10~ 3 < 
Am 2 < 1 x 10 _1 cV 2 is consistent with the observed zenith angle shape at 90% C.L. This result supports 
the evidence for neutrino oscillations given by the analysis of the contained atmospheric neutrino events 
by Super-K. 

We gratefully acknowledge the cooperation of the Kamioka Mining and Smelting Company. The Super- 
Kamiokande experiment has been built and operated from funding by the Japanese Ministry of Education, 
Science, Sports and Culture, and the United States Department of Energy. 
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TABLE I. List of theoretical uncertainties in the flux calculation. 



Error source 


Error (%) 


Chemical composition of the rock 


«l a 


v flux normalization 


±20 a 


Theoretical model dependence 




absolute flux 


±10 a 


bin by bin 


-3.7 to +1.6 b 


spectral index 


±1.4 a 



a Theoretical bin-by-bin uncertainty 
b Theoretical uncorrelated correlated uncertainty 



TABLE II. List of experimental systematic errors in the flux measurement. 



Error source 


Error (%) 


Uncertainty in A9 rec 


< ±l a 


Detection efficiency 


< ±1.2 b 


7m track length cut 


±0.5 C 


Live time 


±0.1 c 


Effective area 


±0.3 C 


PMT gain 


<l c 


Water transparency 


<l c 



Experimental uncorrelated systematic error specific in the most horizontal bin — 0.1<cos6<0 
b Experimental uncorrelated systematic error 
c Bin-by-bin correlated experimental systematic errors 
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FIG. 1. Zenith angle distribution of through-going muons near the horizon observed by Super-K. Filled triangles 
(open circles) indicate events coming from the 180° azimuthal region where the rock overburden is thick (shallow). 
Most of the downward-going (cos 6 > 0) muons denoted by filled triangles are induced by atmospheric neutrinos. 
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FIG. 2. Upward through-going muon flux observed in Super-K as a function of zenith angle. The error bars 
indicate uncorrelated experimental systematic plus statistical errors added in quadrature. The solid histogram 
shows the expected upward through-going muon flux with normalization (a M = —14%) based on the Bartol 
neutrino flux for the null neutrino oscillation case. Also shown as a dotted line is the expected flux assuming the 
best fit parameters at (sin 2 29, Am 2 ) = (0.95, 5.9 x 10~ 3 eV 2 ), q m = +12% for the <-> v T oscillation case. 
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FIG. 3. Allowed region contours at 68% (dotted contour), 90% (thick solid), and 99% (dashed) C.L. obtained 
by the Super-K upward through-going muon analysis on the (sin 2 28, Am 2 ) plane for the <-> v T oscillation 
hypothesis. The star indicates the best fit point at (sin 2 29, Am 2 ) = (0.95,5.9 x 10~ 3 eV 2 ). Also shown is the 
allowed region contour (thin solid) at 90% C.L. by the Super-K contained event analysis. The allowed regions are 
to the right of the contours. 
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